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1. Introduction 


Sustainable energy utilization has become a prominent policy 
target in many countries across the world [1,2]. However, there are 
many concerns over the adoption of various energy resources, such 
as energy-shortage risk, greenhouse gas emission, and renewable 
energy development, posing great challenges for decision mak- 
ers in multi-scale energy management systems (EMSs). Also, the 
management efforts are complicated by climate change due to 
its significant effects on statistics of many energy sources (e.g., 
max, min, mean values and variance of hydropower, solar and 
wind energies). Moreover, such efforts are subject to multiple 
levels of uncertainties that are associated with the forecast of 
changing climatic conditions (e.g., uncertain shifts of precipitation 
patterns), the interactions between energy resources/technologies 
and climatic variations (e.g., uncertain fluctuation of hydropower 
generation in responding to precipitation variations), and the plan- 
ning of EMSs (e.g., uncertain estimation of energy prices). This leads 
to a variety of complexities in relevant decision-making under cli- 
mate change. These complexities may be further multiplied by 
multi-period, multi-facility and multi-objective features as well as 
system dynamics. Thus, system analysis techniques are desired to 
assist in formulating long-term energy management and allocation 
plans, which will be helpful not only for analyzing tradeoffs among 
various socio-economic, energy-related, and environmental objec- 
tives, but also for providing cost-effective adaptation schemes in 
responding to varied impact levels of climate change. 

Previously, a plenty of efforts were undertaken for investigat- 
ing impacts of climate change on multi-scale energy management 
systems and planning the corresponding adaptation capacities 
[3-13,74-78]. For instance, Peramunetilleke [5] assessed impacts 
of climate change on Sri Lanka’s hydroelectric production. They 
believed that power supply from hydropower without back-up 
sources would be highly risky under extreme climatic events. 
Allen and Christensen [14] briefly reviewed the impacts of cli- 
mate change and suggested that a diverse range of issues would 
contribute to this process. Sutherland [15] qualitatively analyzed 
potential impacts of climate change on a number of energy- 
intensive industries in US. Breslow and Sailor [16] investigated 
impacts of climate change on wind speeds and, hence, on wind 
power generation in US. A series of GCM outputs from the Cana- 
dian Climate Center and the Hadley Center were used for providing 
a range of possible variations in seasonal mean wind magnitude. 
Potential impacts of climate change on heating and cooling energy 
demands were investigated through building energy simulation 
models under hourly weather scenarios and were applied to the 
Zurich-Kloten area of Switzerland [17]. Gaterell and McEvoy [10] 
examined uncertain effects of climate change on performance of 
energy insulation measures. Christenson et al. [8] investigated 
impacts of climate change on energy demand of residential build- 
ings through the degree-days method. A procedure to estimate 
heating and cooling degree-days from monthly temperature data 
was proposed, verified, and applied to four representative buildings 
in Switzerland. Mirasgedis et al. [18] adopted a regional climate 
model (i.e., PRECIS) to predict future climatic conditions under dif- 
ferent emissions scenarios. The predicted data were used as inputs 
to a multiple-regression model to examine sensitivity of electricity 
demand in Greece’s interconnected power system under various 
climatic and socio-economic conditions. Pryor and Barthelmie [19] 
qualitatively assessed impacts of climate change on wind energy 
sources, particularly on operation and maintenance of several wind 
farms. Based on a qualitative method, Whitmarsh [20] investigated 
various responses to the impacts of climate change in many energy 
Management systems. Guan [21] attempted to evaluate impacts 
of climate change on building environment through combination 
of weather-data forecasting and building simulation techniques. 


Bassi and Baer [22] employed an integrated modeling approach 
to carry out a country-wide, cross-sector analysis of the interac- 
tions among Ecuador’s energy, social, economic and environmental 
sectors under climate change. Isaac and van Vuuren [23] assessed 
potential energy consumptions for future residential heating and 
air conditioning under climate change. 

Under a series of climate change impact levels, adaptation plans 
of EMSs need to be systematically considered by decision mak- 
ers. Over the past two decades, a large number of studies were 
undertaken on adaptation planning toward climate change in envi- 
ronmental and energy management systems planning [24-50]. 
For example, Huang et al. [51] proposed an inexact-fuzzy multi- 
objective programming model for supporting adaptation planning 
of land resources management in Mackenzie Basin under changing 
climatic conditions. Stakhiv [52] suggested that policy makers and 
water resource managers should be aware of the evolving informa- 
tion on climate change impacts as an activity that was preparatory 
to sound decision making on current water resource management 
actions. Tol et al. [24] made a thorough investigation on studies 
related to climate change adaptation. They suggested that there 
should be four categories of adaptation measures toward climate 
change, including: no adaptation, arbitrary adaptation, observed 
adaptation, and modeled adaptation (optimization). Since climate 
change poses significant challenges for water resource manage- 
ment in Canada, de Loé et al. [53] discussed issues related to 
identification of adaptation measures to water supply systems 
in the near future. Ness et al. [33] examined the adaptation of 
infrastructural institutions under climate change in Norway. Two 
municipalities were used as examples for investigating institu- 
tional responses to floods under climate change. Dessai and Hulme 
[54] presented an assessment framework for identifying adap- 
tation strategies that were robust (i.e., insensitive) to climate 
change uncertainties. The framework was applied to a real-world 
case for supporting water resources management in eastern Eng- 
land. Osbahr et al. [55] explored cross-scale dynamics in coping 
and adjusting responses toward climate change based on qualita- 
tive data from a real-world case in Mozambique. More recently, 
Hoffmann et al. [56] empirically examined the determinants of 
adaptation measures/policies toward climate change. They also 
analyzed the influence of possible climate change effects on the 
scope of adaptation measures through a number of linear regres- 
sion models. Bryan et al. [41] analyzed the adaptation strategies of 
farmers in several countries and analyzed the factors influencing 
their decisions. Deressa et al. [57] identified major strategies used 
by farmers to adapt to climate change in Nile Basin of Ethiopia. They 
also analyzed the factors that might affect their choices of method 
and the barriers to adaptation. 

Particularly, a number of studies were undertaken for exam- 
ining potential adaptation strategies of energy sectors in many 
countries under climate change. For instance, Simeonova [58] com- 
plied and analyzed a variety of policies and measures for adapting 
to climate change in the energy sectors of many central and east- 
ern European countries. Matondo et al. [59] examined the impacts 
of climate change on hydrological regimes, water resources man- 
agement, and hydropower generation. Ruth and Lin [60] explored 
potential impacts of climate change on consumptions of natural 
gas, electricity, and heating oil by residential and commercial sec- 
tors in the state of Maryland of the US. Time series analysis was 
then used to quantify relationships between historical tempera- 
ture and energy demand. A dynamic computer model based on 
those relationships was used for simulating future energy demand 
under a range of energy prices, temperatures, and other driving 
factors. Mansur et al. [61] adopted a national energy model of fuel 
choice to generate climate change adaptation strategies at both 
households and firm scales. They suggested that increasing elec- 
tricity consumption for space cooling would be observed with a 
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reduction in energy consumption for space heating under climate 
change. Overall, energy-related expenditures would likely increase 
in America due to the impacts of climate change on residential 
heating and cooling, resulting in potential welfare damages. Jenk- 
ins et al. [62] attempted to quantify how climate change will have a 
direct effect on heating and cooling energy use in future office envi- 
ronments. The results confirmed the importance of demand-side 
management before assessing the supply-side opportunities under 
extreme climatic events. The study also highlighted the importance 
and possibilities, of adaptation strategies to future climates, and 
the benefits of promoting heating-dominated buildings instead of 
cooling-dominated ones. Weisser et al. [63] used a qualitative anal- 
ysis method to explore the conditions under which nuclear power 
could adapt to climate change over a long term period. 

Previous studies on climate change impact analysis related to 
EMSs are mostly based on qualitative methods, and, thus, could 
merely generate qualitative information. Many studies were con- 
ducted to analyze the impacts of climate change on energy sectors 
as well as to examine the corresponding adaptation strategies 
through integrating many quantitative methods and/or GCM/RCM 
results. However, direct incorporation of GCM/RCM results into 
climate-change impact analysis for multi-scale EMSs was limited to 
one or several climatic sensitive parameters such as temperature, 
precipitation, and humidity, which might be unreasonable in real- 
istic problems. In addition, most of the previous studies focused on 
impacts of climate change on an individual energy sub-sector, pro- 
cess or technology (e.g., residential, commercial and hydropower 
generation sub-sectors). There was a lack of studies that could 
comprehensively reflect and analyze impacts of climate change on 
multiple energy sub-sectors and technologies, and thus analyze the 
impacts from an entire EMS point of view. Moreover, the process 
of impact analysis involved multiple levels of imprecise informa- 
tion, which was largely neglected in the previous studies. In fact, 
the process of climate change impact analysis is subject to a vari- 
ety of uncertainties. For example, in most of the quantitative or 
semi-quantitative methods based on GCM/RCM results, there are 
multiple formats of uncertainties embedded within not only the 
input information, but also the weather-forecasting mechanisms. 
Furthermore, impact analysis could merely provide answers to the 
question of “which sector/process within an energy management 
system would be the most vulnerable under climate change”. Deci- 
sion makers may need to know “if there is any necessity to take 
adaptation actions” and “what the optimal adaptation actions for 
the energy sector under multiple system uncertainties will be”. The 
questions need to be answered through EMSs adaptation planning. 
Likewise, in the process of adaptation planning based on optimiza- 
tion models, relevant parameters and coefficients may be subject 
to multiple formats of uncertainties (e.g., intervals, probabilistic 
and possibilistic distributions). Until now, there were no reported 
studies that could link adaptation planning with climate-change 
impact analysis and could address multi-level uncertainties asso- 
ciated with impact analysis as well as the subsequent adaptation 
planning process. Therefore, an integrated modeling system for 
seamlessly coupling climate change impacts with EMSs adapta- 
tion planning is desired, particularly for a region that is heavily 
dependent on renewable energy resources and is, thus, vulnerable 
to climate change. 


2. Methodology 


2.1. A systematic perspective on adaptation planning under 
climate change 


Almost all energy-related processes and factors would be 
directly and/or indirectly influenced by changing climatic condi- 
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Fig. 1. Climate change impacts and adaptation schemes of EMSs. 


tions (Fig. 1). Among them, end-user demands, renewable energy 
resources and the corresponding utilization facilities would prob- 
ably be the most sensitive and vulnerable ones. For instance, 
availabilities of renewable energy resources (e.g., hydropower and 
wind energy) are highly dependent on the statistics (e.g., min, 
peak and mean values as well as variance) of many meteorolog- 
ical variables that might fluctuate within a certain range due to 
climate change. Such variations of renewable energy availabilities 
would then affect operating statuses of relevant facilities, result- 
ing in changes in their energy outputs. Similarly, in responding 
to variations of climatic factors/conditions, such as temperature, 
solar radiation, and humidity, end-user demands for space heat- 
ing, water heating, and space cooling would correspondingly vary 
over a certain range. Apparently, these impacts would not affect 
an individual energy sub-sector (e.g., energy demand and produc- 
tion) in a separate way. Instead, they interact with each other and, 
thus, would cause an integrated impact on the entire EMS. Such 
an integrated impact over a long-term planning period needs to 
be systematically analyzed by decision makers in order to discern 
the most optimal adaptation schemes for the entire EMS under 
climate change. However, such a task is difficult due to a vari- 
ety of uncertainties associated with changing climatic conditions 
as well as decision makers’ subjective judgments on impacts. For 
example, due to the lack of sufficient information related to the 
evaluation guidelines (e.g., impact levels), fuzziness is always asso- 
ciated with the process of impact identification and integration, 
which is normally caused by vague and imprecise linguistic terms. 
Correspondingly, the evaluation of occurrence chances of climate 
change and their consequences on EMSs is also subject to uncertain- 
ties and is difficult to be handled through conventional methods. 
This implies that planning activities for EMSs would be misled into 
a deviated or even false direction without a comprehensive consid- 
eration over the integrated impacts of climate change. 

Moreover, there is a variety of uncertain information associated 
with decision process, which can only be expressed as interval 
numbers without known distributions. Such uncertain informa- 
tion need to be effectively incorporated into the planning process 
and the obtained decision alternatives, which represents a differ- 
ent level of uncertainties compared with those implied in climate 
change analysis. Therefore, in this research, the integrated climate 
change impacts will be incorporated into the modeling framework 
to facilitate identifying adaptation plans with increased robustness 
under various changing conditions. Specifically, a fuzzy-interval 
inference method (FIIM) and an inexact energy model (IEM) will 
be developed and then integrated into a general modeling system 
(IEMS) as well as applied to the province of Manitoba (IEMS- 
Manitoba, presented in a companion paper). 
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Fig. 2. Membership function of climate change impacts on energy supply (i.e., 
renewable availability factors based on renewable energy resources). Note: 
AF = availability factor. 


2.2. Coupled climate-change impact analysis under uncertainty 


2.2.1. Impacts on energy supply 

Considering an EMS wherein electricity production mainly 
occurs based on renewable energy resources such as hydropower 
and wind energy, over a long-term period, climate change will have 
direct impacts on the availability of these resources due to potential 
variations in many climatic factors and conditions such as precipi- 
tation types, radiation intensities, and wind speeds. Consequently, 
the changes in resource availability would affect operation sta- 
tuses of relevant electricity generation facilities. In this research, 
an availability factor is introduced, which is defined as the amount 
of time that a hydroelectric power station (or wind farm) needs 
to produce electricity over a certain period, divided by the total 
time of this period. Then, five fuzzy sets are assigned to repre- 
sent different levels of impacts of climate change on the availability 
factors of renewable-energy-based power generation facilities, cor- 
responding to “low impact”, “low-to-medium impact”, “medium 
impact”, “medium-to-high impact” and “high impact”. The con- 
struction of membership functions for these five fuzzy sets relies 
on the expertise and experience of experts and stakeholders. The 
involved experts and stakeholders should have in-depth knowl- 
edge about the province and the energy management system, and 
can, thus, provide valuable inputs for quantifying the uncertainties 
[64,65]. Fuzzy logic is then employed for processing such subjective 
opinions. 

In this study, fuzzy sets corresponding to different levels of 
climate change impacts were obtained through questionnaire sur- 
veys. Forty-eight persons were interviewed about their opinions on 
the corresponding relations between climate change impact lev- 
els and renewable energy availability factor decrements. As listed 
in Table 1, 54.10% of the surveyed respondents believed that “low 
impact” should correspond to the statement “the availability fac- 
tor of a renewable-energy-based power generation facility will 
decrease by 5.00% or less”; 58.40% of the respondents selected: “the 
availability factor of a renewable-energy-based power generation 
facility will decrease by approximately 10%” as “low-to-medium 
impact”; 58.40% of the respondents agreed that “the availabil- 
ity factor of a renewable-energy-based power generation facility 
will decrease by around 15.00%” represents “medium impact”; 
60.40% of the respondents matched “the availability factor of a 
renewable-energy-based power generation facility will decrease 
by approximately 20.00%” with “medium-to-high impact”; and 
77.10% of the respondents considered “the availability factor of a 
renewable-energy-based power generation facility will decrease 
25.00% or more” as “high impact”. According to Hwang and Chen 
[66], the membership functions of these fuzzy sets can then be 


Table 1 
Impacts of climate change on energy supply (power generation from renewable 
energy resources). 


The availability factor of renewable 
energy resources will decrease by 
approximately: 


No. of respondents Percentage (%) 


Low impacts 

1% or less 0 0.00 
2% or less 0 0.00 
3% or less 13 27.10 
4% or less 9 18.80 
5% or less 26 54.10 
6% or less 0 0.00 
7% or less 0 0.00 
8% or less 0 0.00 
9% or less 0 0.00 
Total no. of respondents 48 100.00 
Low-to-Medium impacts 

5% 0 0.00 
6% 0 0.00 
7% 0 0.00 
8% 10 20.80 
9% 10 20.80 
10% 28 58.40 
15% 0 0.00 
20% 0 0.00 
25% 0 0.00 
Total no. of respondents 48 100.00 
Medium impacts 

5% 0 0.00 
6% 0 0.00 
7% 0 0.00 
8% 0 0.00 
9% 10 20.80 
10% 10 20.80 
15% 28 58.40 
20% 0 0.00 
25% 0 0.00 
Total no. of respondents 48 100.00 
Medium-to-High impacts 

5% 0 0.00 
6% 0 0.00 
7% 0 0.00 
8% 0 0.00 
9% 0 0.00 
10% 6 12.50 
15% 13 27.10 
20% 29 60.40 
25% 0 0.00 
Total no. of respondents 48 100.00 
High impacts 

5% or greater 0 0.00 
6% or greater 0 0.00 
7% or greater 0 0.00 
8% or greater 0 0.00 
9% or greater 0 0.00 
10% or greater 0 0.00 
15% or greater 0 0.00 
20% or greater 11 22.90 
25% or greater 37 77.10 
Total no. of respondents 48 100.00 


constructed based on above investigation results. The constructed 
fuzzy sets are displayed in Fig. 2. In this figure, “L”, “L-M”, “M”, 
“M-H” and “H” represent low, low-to-medium, medium, medium- 
to-high, and high impacts, respectively. 
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Table 2 

Climate change impacts on end-user energy demands (heating and cooling). 
Energy demands for heating and No. of Percentage 
cooling will increase by approximately: respondents (%) 
Low impacts 
1% or less 23 47.90 
1.5% or less 14 29.20 
2.0% or less 11 22.90 
2.5% or less 0 0.00 
3.0% or less 0 0.00 
3.5% or less 0 0.00 
4.0% or less 0 0.00 
4.5% or less 0 0.00 
5.0% or less 0 0.00 
Total no. of respondents 48 100.00 
Low-to-Medium impacts 
1% 0 0.00 
1.5% 6 12.50 
2.0% 19 39.60 
2.5% 15 31.20 
3.0% 7 14.60 
3.5% 0 0.00 
4.0% 0 0.00 
4.5% 0 0.00 
5.0% 0 0.00 
Total no. of respondents 48 100.00 
Medium impacts 
1% 0 0.00 
1.5% 0 0.00 
2.0% 2 4.20 
2.5% 10 20.80 
3.0% 21 43.80 
3.5% 11 22.90 
4.0% 4 8.30 
4.5% 0 0.00 
5.0% 0 0.00 
Total no. of respondents 48 100.00% 
Medium-to-High impacts 
1% 0 0.00 
1.5% 0 0.00 
2.0% 0 0.00 
2.5% 1 2.00 
3.0% 8 16.70 
3.5% 14 29,20 
4.0% 21 43.80 
4.5% 4 8.30 
5.0% 0 0.00 
Total no. of respondents 48 100.00 


High impacts 

1% or greater 0 0.00 

1.5% or greater 0 0.00 

2.0% or greater 0 0.00 

2.5% or greater 0 0.00 
0 
1 


3.0% or greater 0.00 
3.5% or greater 2.00 
4.0% or greater 3 6.30 
4.5% or greater 21 43.80 
5.0% or greater 23 47.90 
Total no. of respondents 48 100.00 


2.2.2. Impacts on energy demand 

Similar to the availabilities of renewable energy resources, 
energy demands are also sensitive to the variations in climatic 
conditions, especially those for heating and cooling in residential 
and commercial sectors [67,68,73]. Similar to the identification of 
climate change impacts on energy supply, the impacts on end- 
user energy demands are also categorized into five levels (i.e., 
low, low-to-medium, medium, medium-to-high, and high) which 


1.0% 1.5% 2.0% 2.5% 3.0% 3.5% 4.0% 4.5% 5% 


End-user demand 


variation 


Fig. 3. Membership function of the climate change impacts on end-user energy 
demand. 


are presented as five fuzzy sets. A survey was conducted to 
obtain the associated membership functions. The energy demands 
were assumed to increase by 1.00-5.00% compared to the orig- 
inal amounts. The survey results, listed in Table 2 suggest that: 
“energy demands for heating and cooling in the residential and 
commercial sectors will increase by approximately 1.00% or less” 
corresponds to “low impact”; “energy demands for heating and 
cooling in residential and commercial sectors will increase by 
approximately 1.50%” corresponds to “low-to-medium impact”; 
“energy demands for heating and cooling in residential and com- 
mercial sectors will increase by approximately 2.00%” corresponds 
to “medium impact”; “energy demands for heating and cooling 
in residential and commercial sectors will increase by approxi- 
mately 2.50%” corresponds to “medium-to-high impact”; “energy 
demands for heating and cooling in residential and commercial sec- 
tors would increase approximately 3.00% or more” corresponds to 
“high impact”. Fuzzy membership functions associated with these 
five impact levels are established according to the above investiga- 
tion results (Fig. 3). 


2.2.3. Coupled impacts on energy management systems 

In order to comprehensively investigate the impacts of climate 
change on the entire EMS, integrated impact levels are derived 
from combining those impacts on both energy supply and end- 
user demand sub-sectors. However, it is very hard to link climate 
impacts on both sides through conventional quantitative methods. 
Therefore, fuzzy logic could be an effective tool for facilitating the 
quantification of such integrated impacts, leading to the develop- 
ment of a fuzzy-interval inference method (FIIM). 

In this study, the determination of the integrated impact levels 
will be based on a series of fuzzy rules from experts and stakehold- 
ers. The impact levels were set to include six categories, namely, 
“low”, “low-to-medium”, “medium”, “medium-high”, “high” and 
“very-high”. The six levels are expressed as six fuzzy sets. Then, the 
fuzzy logic operator “AND” is used to join factors in the antecedents 
of the rules. Since both environmental and health risks include 
five categories of fuzzy events, there will be a total of 150 rules 
(5x 5x6). Assume that, if a rule achieves the highest frequency 
in the survey, then it will be kept in the rule base for the deter- 
mination of the integrated impact level. It was found that 68.30% 
of the surveyed respondents selected “if both impacts on demand 
and supply are low, then the integrated impact will be low”; 78.50% 
selected “if the impacts on supply side is low and demand side is 
low-to-medium, then the integrated impact level will be low-to- 
medium”; 76.80% chose “if the impact on energy supply is low and 
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Table 3 
Rules for the coupled climate change impact level. 
No. Antecedent Then the integrated impact is 
If impacts on renewable energy availability factor is And impacts on end-user demands is L L-M M M-H H VH 
1 Low Low J 
2 Low Low-to-Medium af 
3 Low Medium af: 
4 Low Medium-to-High af 
5 Low High af 
6 Low-to-Medium Low af 
7 Low-to-Medium Low-to-Medium af: 
8 Low-to-Medium Medium af 
9 Low-to-Medium Medium-to-High gf 
10 Low-to-Medium High NA 
11 Medium Low J 
12 Medium Low-to-Medium af 
13 Medium Medium af 
14 Medium Medium-to-High ri 
15 Medium High J 
16 Medium-to-High Low aj 
17 Medium-to-High Low-to-Medium NA 
18 Medium-to-High Medium VA 
19 Medium-to-High Medium-to-High gf 
20 Medium-to-High High af 
21 High Low af 
22 High Low-to-Medium v 
23 High Medium af’ 
24 High Medium-to-High af 
25 High High af 


Note: L=Low; L-M=Low-to-Medium; M = Medium; M-to-H = Medium-to-High; H=High; VH=Very High. 


the one on energy demand is medium, then the integrated impact 
will be medium”; 83.60% selected “if the impact on energy supply 
is low and the one on energy demand is medium-to-high, then the 
integrated impact will be medium-to-high”; 70.50% selected “if the 
impact on energy supply is low and the one on energy demand is 
high, then the integrated impact will be high” (Table 3). 

Since the integrated impact level can be categorized into “low”, 
“low-to-medium”, “medium”, “medium-high”, “high” and “very- 
high”, the corresponding membership functions of these fuzzy 
events can then be established according to Hwang and Chen [62] 
(Fig. 4). The range of the integrated impact levels is subjectively 
assigned with a range of [0,100]. Thus, they will have single numeri- 
cal site scores after defuzzification. These numerical values have no 
direct relationship with the values of the input impact factors (e.g., 
climate change impact levels of both energy supply and demand 
sides). However, after establishing the fuzzy sets of the integrated 
impact level, a numerical site score can be obtained from Fig. 4 
through the fuzzy “AND” or fuzzy “OR” operations based on the 
two-side impact levels. This will be illustrated in the following case 
study. The site management decisions can then be made based 


0 10 20 30 40 50 60 70 80 90 100 Coupled impact 


Fig. 4. Membership functions of integrated climate change impact levels on energy 
management systems. 


on the calculated site scores that describe the general risk level. 
Table 4 lists the relationship between site scores and suggested 
management actions. 


2.3. Adaptation planning under uncertainty 


The results of integrated climate change impact analysis can 
only provide information related to “how serious will climate 
change affect the EMS?” to decision makers. Decision makers prob- 
ably need to know “how and when to take adaptation actions 
toward climate change?” Therefore, it is desired to couple the 
impact analysis results with adaptation planning efforts. Also, 
uncertainties associated with parameter/coefficients of adaptation 
planning need to be dealt with. Thus, an inexact optimization 
approach will be adopted. In this research, an inexact energy model 
(IEM) will be developed for supporting long-term energy manage- 
ment systems planning under the business as usual (BAU) case, 
as well as adaptation planning under climate change. A number 
of issues concerning cost-effective allocation of energy resources, 
technologies, activities and investments can, thus, be addressed. 
It is expected that the results could provide useful information 
and decision support for energy management systems planning 
within a provincial context. IEM is a large-scale dynamic pro- 
gramming model for energy systems planning under uncertainty 
and complexity. This model is subject to a number of constraints, 
such as end-user demand, economic cost and energy availability. It 


Table 4 
General adaptation strategies under coupled climate change impacts. 


Calculated score Adaptation planning measures 


90-100 Immediate actions need to be taken for mitigating climate 
change impacts on energy management systems 

70-90 Take full actions for mitigating climate change impacts on 
energy management systems 

50-70 Restrict energy export and consumption 

30-50 Take interim control measures 

10-30 Effects of climate change on energy management systems 


should be monitored 
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can be used to describe detailed characteristics of energy-related 
technologies and evaluate economic and environmental effects of 
relevant policies and regulations. Such an energy management 
system can be characterized as a network/chain of energy flows, 
starting from energy-supply options (such as energy production, 
export, and import) and ending in end-use sectors (such as resi- 
dential, commercial, industrial, and transportation sectors). It can 
also be customized into several sub-sectors based on site-specific 
features when it is applied to real-world cases such as in Province 
of Manitoba. 

With consideration of local energy management policies, the 
objective function is to minimize the total costs associated with 
energy services, activities and investments in the study province 
over the entire planning horizon. This objective function is a com- 
bination of costs associated with a variety of processes that move 
energy from production and conversion to consumption. The cost 
of each process is further divided into sub-costs, such as costs 
of energy resource supply, energy process/conversion, technology 
investment, and environmental management. Decision variables 
in IEM represent key discrete points of energy activities. Elec- 
tricity generated from different facilities, production capacities of 
energy generation facilities, volumes of oil imported, and gaso- 
line consumed by private automobiles can be examples of discrete 
variables. The constraints include the relationships among decision 
variables, energy production efficiencies, capital costs, investment 
costs, and energy availabilities. Also, in order to address uncertain- 
ties associated with many factors/processes within an EMS and can 
be expressed as interval numbers, an inexact programming method 
is adopted. The model is based on the interval-parameter linear 
programming (ILP) model, in which the uncertain parameters are 
expressed as intervals without any distributional information that 
is always required in fuzzy and stochastic programming. The ILP 
allows the interval information to be directly communicated into 
the optimization process and resulting solution [69]. According to 
Huang et al. [70-72], an ILP model can be written as follows: 


Minf* = C*X (1a) 
subject to: 

A*X* <B (1b) 
X*>0 (1c) 


where Ate {RE prn, B+ € {RE p1, C € {RE Hn, Xte {R*}@1, 
and R* denotes a set of interval numbers. 
Thus, IEM can be formulated as follows: 


I 
Objective function : Minf* = 5 cE (2a) 


(1) The total cost of energy consumption in the residential sub- 
sector 


T 
ct aa a 

X Cit =a t2 tG tG ts (2b) 

t=1 


(1a) Sub-cost of natural gas consumption in the residential sub- 
sector 


ean 
z= 


T 
1 1 + 1 mn 
t= SP eer cere +X att Xi1,5,t +4116.) (2c) 
t=1 


(1b) Sub-cost of fuel oil consumption in the residential sub- 
sector 


T 
= SOR Oto at +Xi93.¢+Xi2,4,0) (2d) 
t1 


(1c) Sub-cost of liquefied petroleum gas (LPG) consumption in 
the residential sub-sector 


T 


+ _ + + + + 
ci3 = ` RE 3 (XT 3,1,t + X133, +X13,4,t + X1,3,5,t + X136, (2€) 
t=1 


(1d) Sub-cost of heating oil consumption in the residential sub- 
sector 


Cha = Gas +Xfaar) (2f) 


(1e) Sub-cost of renewable energy resources consumption in the 
residential sub-sector 


T 
c= She iG oa + Xs oetXt sac +%1.5,6.0) (28) 
t1 


(1f) Sub-cost of electricity consumption in the residential sub- 
sector 


T 
oa fi + fi = fi cai 1 an 

Cigs N RE eXe t XT6,3,t + X1,6,4,t + X1,6,5,t +X1,6,6,.(2h) 
t=1 


(2) The total cost of energy consumption in the commercial sub- 
sector 


T 

i m f 
` Gt = 1 +2 + CO3 + O4 t ast 6 (2i) 
t1 


(2a) Sub-cost of natural gas consumption in the commercial sub- 
sector 


T 
+ + + : 
GiS WS Oan + X313, + X31,6,0) (2j) 

t=1 


(2b) Sub-cost of fuel oil consumption in the commercial sub- 
sector 


T 
Za = Y R221 + X2,3,0 tX at tX) (2k) 
t=1 


(2c) Sub-cost of LPG consumption in the commercial sub-sector 


T 

= aa + fi 1 fi = 

Coz X Ra Oat EX33,3,t + X3,3,5,t tX2,3,6,0) (21) 
t=1 


(2d) Sub-cost of heating oil consumption in the commercial sub- 
sector 


T 


+ + 
Bas X RG aant +X3 4.3.1) 
t1 


(2m) 


(2e) Sub-cost of renewable energy resources consumption in 
commercial sub-sector 


T 


+ + + 
O5= X ORs Osis + X3'5,3,t +X.5,6,1) (2n) 
t=1 


(2f) Sub-cost of electricity consumption in commercial sub- 
sector 


T 
+ + 
CS ` RE 6 t(X36,1,t t XZ6,2,t + X26,3,t + X26,4,t + X26,5,t 


+X36,6,t) (20) 
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(3) The total cost of energy consumption in the transportation 
sub-sector 


T 
+ + 

` Gt = 031 HG2 + CG3 

t=1 


(3a) Sub-cost of natural gas consumption in the transportation 
sub-sector 


E 1 


3.47 (2p) 


+ 
C35 +036 +037 


T 


+ + + 
Ce ` Rt XS tX31,2,t +X 132+ X31, 
t=1 


+X3 45.4 +X3 160) (2q) 


(3b) Sub-cost of fuel oil consumption in the transportation sub- 
sector 


T 


+ + 
i i X REG 9 (XZ 05,6 + X32,6,t + X3.2,7,0) 
t=1 


(2r) 


(3c) Sub-cost of LPG consumption in the transportation sub- 
sector 


T 


Hee a + 
er ` Res (X33. 16 + X332, + X33311 
t=1 


xt ) 


3,3,5,t (2s) 


(3d) Sub-cost of diesel consumption in the transportation sub- 
sector 


T 


Foai + 
Ca X RE 4 (X34,1,t +X3,4,2,t + X3,4,3,t +X3,4,5,t +X3,4,6,t) (2t) 
t=1 


(3e) Sub-cost of jet fuel consumption in the transportation sub- 
sector 


T 
CoS ` RE 5 X354, 
t=1 


(3f) Sub-cost of gasoline consumption in the transportation sub- 
sector 


(2u) 


T 


— + 
03,6 = X RE 6,t(X3,6,1,t + X3,6,2,t + X3,6,3,t + X3,6,5,t + %36,6, (2V) 
t=1 


(3g) Sub-cost of electricity consumption in the transportation 
sub-sector 


+ 
x X375, 7 


3,7,4,t T X3, 


3,7,6, (2w) 


T 

A + 

G7 ` RE 9 (XS 71 +X37,3,t7 
t=1 


(4) The total cost of energy consumption in the industrial sub- 
sector 


T 

at + 
SGi = C41 + C42 + C43 1 
t=1 


(4a) Sub-cost of natural gas consumption in the industrial sub- 
sector 


M T 
CaS ` 5 Rtgs Xam. 
t=1 


m 


Chat C45 (2x) 


(2y) 


(4b) Sub-cost of fuel oil consumption in the industrial sub-sector 


M 


T 
2 = ` ` RIZ 2,tX4.2,m,t 
t=1 


m 


(2z) 
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(4c) Sub-cost of LPG consumption in the industrial sub-sector 


M 


GIS ` 


m 


T 
X RIZ 3,tX4.3,m,t 
t= 


(2aa) 


(4d) Sub-cost of coal consumption in the industrial sub-sector 


T 
Cha = ` 
t=1 


X, 


£ 


X, 


HX. 


Rega (X34111 


4,4,3,t + 


4,4,5,t | 


4,4,6,t) 


(2ab) 


(4e) Sub-cost of renewable energy resources consumption in the 


industrial 


M 
5 = ` 


m 


sub-sector 


T 
X Rtgs tX4.5,m,t 
t=1 


(2ac) 


(4f) Sub-cost of electricity consumption in the industrial sub- 


sector 
M 


C46 = X 


m 


T 
X REE emi 
t=1 


(5) The total cost of power generation facilities 


| (AFOM} į + AVOM;;_Y, 


n, 


(2ad) 


Fr) + CAPS INVE e Ze] 


n,t “n,t 
t=1 n=1 t=1 
(2ae) 
subject to: 
(1) Balance of energy carriers 
(1a) Balance of natural gas 
2 5 7 5 
Se t X Xaime + SOG ia 
i=1 m=1 m=1 m=1 
( Eley, ) y (2af) 
+{ aoe | SMi ve a 
Conv; , i 
(1b) Balance of fuel oil 
ate fi fi fi i $ 1 1 
Xia + X1,2,3,t + X1,2,4,t + X2,2,1,t + %2,2,3,t + X2,2,4,¢ + X2,2,6,t 
6 
1 i fi $ I ob + £ 
X325, + X3,2,6,t +%3,2,7,.0 7 DTT, < N3, IMP Ve 
m=1 


(2ag) 
(1c) Balance of LPG 
p= i £ i | E T i 
X1,3,1,t t X1,3,3,t + X1,3,4,t +¥7,3,5,¢ + %1,3,6,t + X2,3,1,¢ + %2,3,3,t 
J | L £ J ii | 
+X9°3.5,¢ + X2,3,6,t + %3,3,1,t t %3,3,2,t t %3,3,3,t t %3,3,5,t 
M T 
+ X Xia mt S15 IMP3 p yt (2ah) 
m t=1 
(1d) Balance of heating oil 
+ + a . 
Xa t X1,4,4,t + X2,4,1,¢ + X2,4,3,¢ < M4, MP3 Ve (2ai) 
(1e) Balance of renewable energy resources 
£ li E p Ji 
X1,5,1,t t X1,5,2,t t X1,5,4,t + %1,5,6,t + %2,5,1,t + %2,5,3,t 42.5.6. 
M T 
+X X RiEs X$s mı < DREN, Vt (2aj) 
t=1 


m 
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(1f) Balance of electricity 


N 
+ EleFO; + EleCOAL; +EleREN; > X` (Y#,), 


EleNGA; Vt (2ak) 
n=1 
(1g) Balance of diesel 
£ | | | E 
X3,3,1,t t %3,3,2,t t %3,3,3,t t %3,3,5,t + %3,3,6,t 
Ast < në, IMPĖ,, Vt (2al) 
Conv, T t >t 
(1h) Balance of jet fuel 
+ 4 
X35 4t < MPE Ve (2am) 
(1i) Balance of gasoline 
X3 6,1,t +X3.6,2, + X3.6,3,¢ +X36,5,t +%3.6,6.c © N7 MP7 Vt 
(2an) 
(1j) Balance of coal 
3t 
Xpatet+Xaa3e+%Xa45,¢+%aaee4 Conve = epee Vt 
t 
(2ao) 
(2) End-user energy demand constraints 
(2a) Space heating demands 
Ne Conf, Xia > DSH;,, Vi,ti=1and2 (2ap) 
(2b) Space cooling demands 
X Cont Xii > DSC, Vi, t, i=1and2 (2aq) 
(2c) Water heating demands 
S “cons X53 2 DWH Vi,t,i=1and2 (2ar) 
(2d) Lighting demands 
Seo $ qe Z DLG}, Vitis 1and2 (2as) 
(2e) Appliance energy demands 
Leonti $ s2 DAPR, Vi, t,i=1and2 (2at) 
(2f) Commercial ventilation demand 
Eont Bar 2 DVP, Viti=2 (2au) 
(2g) Light duty vehicle energy demand 
Sows Xia, 2 DLD}, VYi=3,t,j=1,3,4,6and7 (2av) 


(2h) Freight truck energy demand 


4 ; . 
X Conve, Xe ot > DFR; Vi=3,t,j=1,3,4,6and7 


(2i) Bus energy demand 


+ : : 
X conve, Xaa > DBS; Vi=3,t,j=1,3,4,6and7 


(2j) Air plane demand 


S “convi Xj 4,2 DARE, Vi=3,tj=5 


(2k) Recreational energy demand 
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(21) Railway energy demand 


Lows XË o. = Del 


(2m) Equivalent lubricant energy demand 


X Conv Xia > DLB, Vi = 3, t, j = 2 for lubricant 


mi Total electricity demand 
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(20) Agricultural energy demand 
J 
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(2p) Forest energy demand 
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(2q) Mining energy demand 


J 
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(2r) Construction energy demand 


Seow X = DCONŻ„, Vi=4,t,j=1,2,3,5,6 
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(2s) Pulp and paper energy demand 


EE ; ; 
S convi Xi 5 > DPP; Vi = 4,t,j=1,2,...,6 
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(2aw) 


(2ax) 


(2ay) 


(2az) 


Vi = 3,t,j = 1, 2, 4, 6, 7 for railway 


(2ba) 


(2bb) 


(2bc) 


(2bd) 


(2be) 


(2bf) 


(2bg) 


(2bh) 
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(2t) Manufacture energy demand 


Lows X5 e2 DMAN},, Yi=4,tj=1,2,...,6  (2bi) 
(2u) Electricity export targets 
EELE¥ > TARGET_EX;, Vt (2bj) 
(3) Electricity generation facility capacity constraint 
Yat < EXCAP;, e + CAP; oe (2bk) 
LowerY7,<Y;,<UpperY;,, n,t (2bl) 
(4) GHG emission constraints 
(4a) Da ij, ae < (EmCO3)F (2bm) 
i=1 j=1 =1 
M 
40) SSS Emre tie) time < EmCHa) (2bn) 
i=1 j=1 m=1 
M 
(4c) Dae ist) Xi m, < (EMN20); (2bo) 


i=1 jel m=1 


where t= 1,2,...,5; c11 =the total cost of natural gas (NGA) con- 
sumption in residential sub-sector; c,,2 =the total cost of fuel oil 
(FO) consumption in residential sub-sector; c13 =the total cost of 
liquefied petroleum gases (LPG) consumption in residential sub- 
sector; c,,4=the total cost of heating oil (HO) consumption in 
residential sub-sector; c1, =the total cost of renewable energy 
(RE) resources consumption in residential sub-sector; c1 =the 
total cost of electricity (ELE) consumption in residential sub- 
sector; C21 =the total cost of NGA consumption in commercial 
sub-sector; cz, =the total cost of FO consumption in commercial 
sub-sector; c23 =the total cost of LPG consumption in commercial 
sub-sector; C24 =the total cost of HO consumption in commercial 
sub-sector; C25 =the total cost of RE consumption in commercial 
sub-sector; C2 =the total cost of ELE consumption in commercial 
sub-sector; c31 =the total cost of NGA consumption in trans- 
portation sub-sector; c32 =the total cost of FO consumption in 
transportation sub-sector; c33=the total cost of LPG consump- 
tion in transportation sub-sector; c34=the total cost of diesel 
(DSL) consumption in transportation sub-sector; c35=the total 
cost of jet fuel (JF) consumption in transportation sub-sector; 
c3 6 = the total cost of gasoline (GSL) consumption in transportation 
sub-sector; c37=the total cost of ELE consumption in trans- 
portation sub-sector; c4, =the total cost of NGA consumption 
in industrial sub-sector; c42 =the total cost of FO consumption 
in industrial sub-sector; c43 =the total cost of LPG consumption 
in industrial sector; c44=the total cost of coal (CL) consump- 
tion in industrial sector; c45=total cost of RE consumption in 
industrial sector; c4g=total cost of ELE consumption in indus- 
trial sector; n= electricity generation facilities (n = 1 =hydropower; 
n=2=coal; n=3=NGA; n=3= fossil fuel); t= coefficient of loss; 
Not = = line loss of electricity transportation. 

Obviously, the above IEM model can effectively reflect: (i) inter- 
actions between a number of sub-sectors, technologies and energy 
resources among various energy management systems, (ii) uncer- 
tain information associated with input parameters that can be 
expressed as intervals, and (iii) the provincial baseline of energy 
management over the planning horizon. Thus, a series of decision 
related to “how and when to take adaptation actions toward climate 
change” can be obtained. 


3. Method of solution 


The method of solution is based on (a) the results of FIIM for 
climate change impact analysis, and (b) the interactive solution 
algorithm of ILP. Such an interactive solution algorithm has been 
developed to solve model (1) through analyzing the detailed inter- 
relationships between parameters and variables and between the 
objective function and constraints. According to the algorithms, 
the solution for interval linear programming (ILP) model can be 
obtained through a two-step method, where a submodel corre- 
sponding to f~ (when the objective function is to be minimized) 
is first formulated and solved, and then the relevant submodel cor- 
responding to f* can be formulated based on the solution of the first 
submodel. In detail, the first submodel can be formulated as follows 
(assume that b> > 0, and f* > 0): 


ky n 
=s + 
Minf~ = ya x; + 5 C %; (3a1) 
j=l j=ky+1 
ai to: 
n 
Sr Sign(a Xj + 5 lay |" Sign(a; x} < bý, Vi (3a2) 
j=ky+1 
xF>0, Vj (3a3) 


J 
k, are interval variables with positive coef- 
ficients in the objective function and x, j=ki+1,k, +2,...,n are 
interval variables with negative coefficients in the objective func- 
tion. Thus, solutions of XT opt (j=1,2,...,k1) and xP opt G=k, +1, 
kı +2,...,n) can be obtained through Salvin submodel (3a). Then 
the submodel corresponding to f* can be formulated as follows 
(assume that b> > 0, and f* > 0): 


where xr, j=1,2,..., 


ky n 
+ ++ + 
Minf De x + 5 CX; (3b1) 
j=l j=ky+1 
subject to: 
ky 
X laji Sign(a; ) xy + 5 lajl" Sign(a; )x; <b, Vi (3b2) 
j=! j=ky4+1 
xe >0, Vj (3b3) 
x; =F opt? j=1,2,...,ky (3b4) 
x, =F opt? j=k +1, ko +2,...,n (3b5) 


Solutions of xt one (j=1,2,...,k1) and XT opt G=k, +1,k; +2,...,.n) 
can be obtained through solving submodel | (3b). Thus, we can have 
the final solution of foo, = [fope»fopr] and x XE opt = PF opt? xy apell 

Thus, the vite ated modeling system (IEMS) can be developed 
through linking the fuzzy-interval inference method (FIIM) and the 
inexact energy model (IEM). The general solution procedure for the 


proposed IEMS is presented in Fig. 5 and described as follows: 


- Step 1: Identify climate change impacts on energy supply based 
on the obtained fuzzy sets; 

- Step 2: Identify climate change impacts on end-user demand 
based on the obtained fuzzy sets; 

- Step 3: Identify multi-layer, multi-sector and multi-period 
impacts of climate change on the EMS through the integration 
method; 

- Step 4: Based on interval analysis and fuzzy sets theory, identify 
the integrated impact levels of climate change on the EMS; 
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Fig. 5. Solution procedure. 


- Step 5: Based on the integrated impact results and the values in 
Table 4, judge whether adaptation actions would be needed for 
the EMS under climate change; 

- Step 6: If no adaptation actions would be needed, formulate the 

inexact model (IEM) under the BAU case; 

Step 7: Transform the IEM model into two submodels, where the 

lower bound of f= is first desired since the objective is to minimize 

f5; 

- Step 8: Formulate the f- submodel, including the f~ objective 

function and the related constraints; 

Step 9: Solve the f submodel according to conventional LP (linear 

programming) solution methods, and obtain the solution for font 

and x7 opt 

- Step 10: Formulate the f* submodel, including the objective func- 

tion and constraints corresponding to f*, as well as the second set 

of constraints regulating the bounds of decision variables; 

Step 11: Solve the f* submodel according to conventional LP (lin- 


ear programming) methods and obtain the solution for opt and 
+ 


X ot! 

Step 12: The optimal results for the IEMS model can, thus, be 

obtained as fópt = [opt fopt] and X opt = PF opt? xy spil: 

Step 13: If adaptation actions need to be undertaken in respond- 

ing to the integrated climate change impacts, formulate the 

corresponding model based on modeling structure formulated 

under the BAU case and the impact analysis results (i.e., the vari- 

ations on energy supplies and demands) obtained from Steps 1 to 

5; 

Step 14: Repeat Steps 7 to 12; 

- Step 15: Make comparisons between results under the BAU case 
and climate change. 


4. Conclusions 


In this research, a large-scale integrated modeling system (IMS) 
was developed for supporting adaptation planning under impacts 


of climate change from an entire EMS point of view. A number 
of methodologies were seamlessly integrated into IMS, includ- 
ing the fuzzy-interval inference method (FIIM), inexact energy 
model (IEM), and uncertainty analysis. Compared to conventional 
studies in the similar areas, such a system could (i) incorporate 
multiple technologies, energy resources, and sub-sectors, and cli- 
mate change impacts within a general modeling framework, (ii) 
address interactions of climate change impacts on multiple energy 
sub-sectors and resources within an EMS, (iii) identify optimal 
adaptation strategies of an EMS to climate change impact through 
a two-step procedure, (iv) deal with multiple levels of uncertain 
information associated with processes of climate change impact 
analysis and adaptation planning, and (v) seamlessly combine 
climate change impact analysis results with inexact adaptation 
planning. Thus, decision makers could have a comprehensive view 
on the EMS as well as the corresponding adaptation schemes under 
climate change, greatly improving the robustness and complete- 
ness of the decision-making processes. The developed method can 
be applied to long-term adaptation planning in energy manage- 
ment systems, particularly those that are highly dependent on 
renewable energy sources. The solutions can generate decision 
alternatives and thus help decision makers identify desired policies 
under climate change. In the companion paper, application of the 
developed approach will be conducted in the Province of Manitoba, 
Canada. 
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